DiFX2: A more flexible, efficient, robust and powerful software correlator

Software correlation, where a correlation algorithm written in a high-level language such as C++ is run on commodity computer hardware, has become increasingly attractive for small to medium sized and/or bandwidth constrained radio interferometers. In particular, many long baseline arrays (which typically have fewer than 20 elements and are restricted in observing bandwidth by costly recording hardware and media) have utilized software correlators for rapid, cost-effective correlator upgrades to allow compatibility with new, wider bandwidth recording systems and improve correlator flexibility. The DiFX correlator, made publicly available in 2007, has been a popular choice in such upgrades and is now used for production correlation by a number of observatories and research groups worldwide. Here we describe the evolution in the capabilities of the DiFX correlator over the past three years, including a number of new capabilities, substantial performance improvements, and a large amount of supporting infrastructure to ease use of the code. New capabilities include the ability to correlate a large number of phase centers in a single correlation pass, the extraction of phase calibration tones, correlation of disparate but overlapping sub-bands, the production of rapidly sampled filterbank and kurtosis data at minimal cost, and many more. The latest version of the code is at least 15% faster than the original, and in certain situations many times this value. Finally, we also present detailed test results validating the correctness of the new code.Comment: 28 pages, 9 figures, accepted for publication in PAS

The Bias and Uncertainty of Redundant and Sky-based Calibration Under Realistic Sky and Telescope Conditions

The advent of a new generation of low frequency interferometers has opened a direct window into the Epoch of Reionisation (EoR). However, key to a detection of the faint 21-cm signal, and reaching the sensitivity limits of these arrays, is a detailed understanding of the instruments and their calibration. In this work we use simulations to investigate the bias and uncertainty of redundancy based calibration. Specifically, we study the influence of the flux distribution of the radio sky and the impact of antenna position offsets on the complex calibration solutions. We find that the position offsets introduce a bias into the phase component of the calibration solutions. This phase bias increases with the distance between bright radio sources and the pointing center, and with the flux density of these sources. This is potentially problematic for redundant calibration on MWA observations of EoR fields 1 and 2. EoR field 0, however, lacks such sources. We also compared the simulations with theoretical estimates for the bias and uncertainty in sky model based calibration on incomplete sky models for the redundant antenna tiles in the MWA. Our results indicate that redundant calibration outperforms sky based calibration due to the high positional precision of the MWA antenna tiles

A microlensing measurement of the size of the broad emission line region in the lensed QSO 2237+0305

We present spatially resolved spectroscopic images of the gravitationally lensed QSO 2237+0305 taken with the GMOS Integral Field Unit (IFU) on the Gemini North telescope. These observations have the best spatial resolution of any IFU observations of this object to date and include the redshifted CIII] and MgII QSO broad lines. Unlike Mediavilla et al. 1998, we find no evidence for an arc of resolved broad line emission in either the CIII] or MgII lines. We calculate the image flux ratios of both the integrated emission lines and the surrounding continua. The flux ratios of the CIII] and MgII emission lines are consistent with each other but differ substantially from their corresponding continuum flux ratios and the radio/mid-IR flux ratios previously published. We argue that the broad emission line region must be microlensed and the CIII] and MgII emission regions must be approximately the same size and co-located along the line-of-sight. Assuming a simple model for the broad emission line region and the continuum region, we show the size of the CIII] / MgII broad line region is ~0.06 h_{70}^{1/2} pc and the continuum region is <= 0.02 h_{70}^{1/2} pc.Comment: 7 pages, 7 figures. Accepted for publication in MNRA

A VLBA search for binary black holes in active galactic nuclei with double-peaked optical emission line spectra

We have examined a subset of 11 active galactic nuclei (AGN) drawn from a sample of 87 objects that possess double-peaked optical emission line spectra, as put forward by Wang et al. (2009a) and are detectable in the FIRST survey at radio wavelengths. The double-peaked nature of the optical emission line spectra has been suggested as evidence for the existence of binary black holes in these AGN, although this interpretation is controversial. We make a simple suggestion, that direct evidence of binary black holes in these objects could be searched for in the form of dual sources of compact radio emission associated with the AGN. To explore this idea, we have used the Very Long Baseline Array to observe these 11 objects from the Wang et al. (2009a) sample. Of the 11 objects, we detect compact radio emission from two, SDSS J151709+335324 and SDSS J160024+264035. Both objects show single components of compact radio emission. The morphology of SDSS J151709+335324 is consistent with a recent comprehensive multi-wavelength study of this object by Rosario et al. (2010). Assuming that the entire sample consists of binary black holes, we would expect of order one double radio core to be detected, based on radio wavelength detection rates from FIRST and VLBI surveys. We have not detected any double cores, thus this work does not substantially support the idea that AGN with double-peaked optical emission lines contain binary black holes. However, the study of larger samples should be undertaken to provide a more secure statistical result, given the estimated detection rates.Comment: 14 pages, 3 figures. To appear in A

A GPU based Transient Dedisersion Search Engine for CRAFT

The Commensal Realtime ASKAP Fast Transient Survey (CRAFT[5]) will search the ASKAP data stream for fast (< 5sec) transient events, associated with the most extreme conditions in the Universe. The CRAFT search will run in parallel with all normal observing on ASKAP, giving a tremendous advance in the transient parameter space which can be searched. ASKAP offers high sensitivity, high resolution, and continous observations of a significant portion of the sky. However, to search across the wide field of view and at the data rates which are being provided by ASKAP is extremely challenging. Nevertheless, as reported here, we are on track to achieve the goals as laid out in the Survey Science Proposal. Using GPUs as a simple highly parallel compute-engine we can monitor the full field of view with a 5(sigma) sensitivity of ~Jy for a millisecond event covering the astronomically significant range of DMs. After that trigger detection we can download the beamformer data-buffer and image the sky at full sensitivity and spatial resolution with an arbitary frequency and time resolution

Astrophysical Supercomputing with GPUs: Critical Decisions for Early Adopters

General purpose computing on graphics processing units (GPGPU) is dramatically changing the landscape of high performance computing in astronomy. In this paper, we identify and investigate several key decision areas, with a goal of simplyfing the early adoption of GPGPU in astronomy. We consider the merits of OpenCL as an open standard in order to reduce risks associated with coding in a native, vendor-specific programming environment, and present a GPU programming philosophy based on using brute force solutions. We assert that effective use of new GPU-based supercomputing facilities will require a change in approach from astronomers. This will likely include improved programming training, an increased need for software development best-practice through the use of profiling and related optimisation tools, and a greater reliance on third-party code libraries. As with any new technology, those willing to take the risks, and make the investment of time and effort to become early adopters of GPGPU in astronomy, stand to reap great benefits.Comment: 13 pages, 5 figures, accepted for publication in PAS

Differential microlensing measurements of quasar broad-line kinematics in Q2237+0305

The detailed workings of the central engines of powerful quasars remain a mystery. This is primarily due to the fact that, at their cosmological distances, the inner regions of these quasars are spatially unresolvable. Reverberation mapping is now beginning to unlock the physics of the Broad Emission Line Region (BELR) in nearby, low-luminosity quasars, however it is still unknown whether this gas is dominated by virial motion, by outflows, or infall. The challenge is greater for more distant, powerful sources due to the very long response time of the BELR to changes in the continuum. We present a new technique for probing the kinematic properties of the BELR and accretion disk of high-z quasars using differential microlensing, and show how substantial information can be gained through a single observation of a strongly-lensed quasar using integral field spectroscopy. We apply this technique to GMOS IFU observations of the multiply-imaged quasar Q2237+0305, and find that the observed microlensing signature in the CIII] broad emission line favours gravitationally-dominated dynamics over an accelerating outflow

Spectroscopy with the Engineering Development Array: cold H+^{+} at 63 MHz towards the Galactic Centre

The Engineering Development Array (EDA) is a single test station for Square Kilometre Array (SKA) precursor technology. We have used the EDA to detect low-frequency radio recombination lines (RRLs) from the Galactic Centre region. Low-frequency RRLs are an area of interest for future low-frequency SKA work as these lines provide important information on the physical properties of the cold neutral medium. In this project we investigate the EDA, its bandpass and the radio frequency interference environment for low-frequency spectroscopy. We present line spectra from 30 to 325 MHz for the Galactic Centre region. The decrease in sensitivity for the EDA at the low end of the receiver prevents carbon and hydrogen RRLs to be detected below 40 and 60 MHz respectively. RFI strongly affects frequencies in the range 276-292, 234-270, 131-138, 95-102 and below 33 MHz. Cn$\alpha$ RRLs were detected in absorption for quantum levels n = 378 to 550 (39-121 MHz) and in emission for n = 272 to 306 (228-325 MHz). Cn$\beta$ lines were detected in absorption for n = 387 to 696 (39-225 MHz). Hn$\alpha$ RRLs were detected in emission for n = 272 to 480 (59-325 MHz). Hn$\beta$ lines were detected for n = 387 to 453 (141-225 MHz). The stacked Hn$\alpha$ detection at 63 MHz is the lowest frequency detection made for hydrogen RRLs and shows that a cold (partially) ionized medium exists along the line of sight to the Galactic Centre region. The size and velocity of this cold H$^{+}$ gas indicates that it is likely associated with the nearby Riegel-Crutcher cloud.Comment: 18 pages, 6 figures and 5 table